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Abstract 

Aims. Up to now, the study of the star formation rate in galaxies has been mainly based on the Ha emission-line luminosity. However, 
this standard calibration cannot be applied at all redshifts given one instrumental setup. Surveys based on optical spectroscopy do not 
observe the Ha emission line at redshifts higher than z ~ 0.5. Our goal is to study existing star formation rate calibrations and to 
provide new ones, still based on emission-line luminosities, which can be applied for various instrumental setups. 
Methods. We use the SDSS public data release DR4, which gives star formation rates and emission-line luminosities of more than 
100000 star-forming galaxies observed at low redshifts. We take advantage of this statistically significant sample in order to study the 
relations, based on these data, between the star formation rate and the luminosities of some well-chosen emission lines. We correct the 
emission-line measurements for dust attenuation using the same attenuation curve as the one used to derive the star formation rates. 
Results. We confirm that the best results are obtained when relating star formation rates to the Ha emission line luminosity, itself 
corrected for dust attenuation. This calibration has an uncertainty of 0.17 dex. We show that one has to check carefully the method 
used to derive the dust attenuation and to use the adequate calibration: in some cases, the standard scaling law has to be replaced 
by a more general power law. When data is corrected for dust attenuation but the Ha emission line not observed, the use of the Hj8 
emission line, if possible, has to be preferred to the [Oil] A 3727 emission line. In the case of uncorrected data, the correction for dust 
attenuation can be assumed as a constant mean value but we show that such method leads to poor results, in terms of dispersion and 
residual slope. Self-consistent corrections such as previous studies based on the absolute magnitude give better results in terms of 
dispersion but still suffer from systematic shifts, and/or residual slopes. The best results with data not corrected for dust attenuation 
are obtained when using the observed [Oll]A3727 and H/3 emission lines together. This calibration has an uncertainty of 0.23 dex. 

Key words. Galaxies: fundamental parameters - Galaxies: statistics - ISM: HII regions - ISM: dust, extinction 



. 1. Introduction 

^ Massive ongoing spectroscopic surveys of the high-redshift 
. Universe allow us to draw a new picture of how the galaxies 
■ evolve through cosmic times. Many physical properties are used 
" as tracers of the evolutionary stage and the past history of one 
galaxy (e.g. stellar mass, metallicity, amount of dust, star forma- 
tion rate}/IlianksJo_ex works on stellar population mod- 
els (iBruzual & Charlotll2003t (pioc & Rocca-Volmerange 1997h 
and photo-ionization models (Chariot & Longhetti 2001), most 
of these parameters can now be well recovered from rest- frame 
optical spectroscopic observations, using all information avail- 
able from the stellar continuum and absorption- or emission-line 
measurements. 

Unfortunately, these models often require a better signal- 
to-noise ratio, wavelength coverage, and/or spectral resolution 
than the ones available on recent deep surveys (e.g. VVDS: 
[Le Fevre et alJ 120051: zCOSMOS: Lily et al. 2007). Thus, a 
simpler approach is to use instead emission-line diagnos- 
tics that would be calibrated on high qua lity samples, taken 
from local Universe surveys (e.g. 2dFGRS: Colless et"aD 120011: 
SDSS: Yorketal. 2000). This approach has already been ex- 
tensively used with older data to recover metallicities with 
oxygen-to-hydrogen (lMcGaughlll99ll: iKewlev & Dopitall2QQ2h 



or ni trogen-to-hydrogen (iPettini & Pagell [20041: Zee et al. 
Il998h line ratios, dust content wit h the Ha/Hj3 B aimer decre- 
ment (ISeatonll 19791 ICalzettill2QQ"l . or star formation rates with 
the Ha or [Qui A 3727 emission-line luminosities (iKennicuttI 
ll998l:lKewlevetal.ll2004 . 

Recently, the star formation rates of > 100000 SDSS 
star-form ing galaxies have b e en e stimated using the CLOl 
method (^Chariot & Longhettil I2QQ1 ) which takes into ac- 
count all in formation available in emission-line measurements 
(iBrinchmann et al.l |2QQ4|) . Then, the emission-line measure- 
ments of the same galaxies have also been made publicly avail- 
able. We can thus use this large amount of data to derive new 
good quality calibrations of the star formation rates (hereafter 
SFR) from emission-line luminosities. 

The standard calibration of the SFR versus the Ha lumi- 
nosity is based on a very simple scaling relation, driven by the 
amount of Ha recombination photons being directly propor- 
tional to the intensity of the ionizing source, i.e. the amou nt of 
young, hot stars, which the SFR is an indicator of (KennicuttI 
Il998|) . The only drawback to this simple assumption is that the 
Ha luminosity has to be corrected for dust attenuation. 

As shown by Kewlev et"^ (|2004j) , the calibration of the 
SFR versus the [Oil] A 3727 luminosity is not only more affected 
by dust attenuation, but also depends on metallicity, since in- 
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terstellar medium with higher gas-phase metallicities produces 
less [On] A 3727 photons because of increased oxygen cooling. 
Starting from this result, they have derived new calibrations of 
the SFR versus [Oil] A 3727 luminosity which take into account 
the metallicity. The idea was to calibrate the [Oii]A3727/Ha 
emission-line ratio versus the metallicity, and then to use it to 
recover the SFR from the Ha luminosity, itself derived from 
the [On] A 3727 luminosity. The dependence on the metallicity 
of the [On] A3727/Ha emission-line ratio, or more generally the 
[Oil] A 3 727 based S FR ca libra tions, has been later confirmed 
by Mo uhcine et al.l (l2005l) and iBicker & Fritze-v. Alvensleb"enl 
(l2Q05i) . 

Recently iMoustakas et al.1 (l2QQ6h and IWeiner et all (l2QQ7l) 
have derived a new calibration of the SFR of galaxies versus 
[On] A 3727 luminosity. They take care of the effect of the metal- 
licity with a slightly different approach, by parametrizing their 
calibration in terms of the B-hand luminosity, which can be 
used as a rough met allicity indicato r because of the metallicity- 
luminosity relation (iLequeux et al.|[l97 9: Skillman et al. 1989; 
iRicher & McCall 1995; Lam areille et al.i i2004i: iTremonti et al.i 
2004; Lamareille et al. 2006). 

We point out that these previous studies base their new cal- 
ibrations on either a direct estimate or a secondary indicator 
(e.g. rest-frame luminosity) of the metallicity, which in both 
cases suffers from intrinsic uncertainties and degeneracies (see 
[Kewlev & ElHson 2008 for a detailed discussion on the uncer- 
tainties of metallicity calibrations, or lLamareille et al.ll2004l for 
a discussion on the dispersion of the luminosity-metallicity rela- 
tion). Moreover, some of these previous works also assume that 
the dust attenuation is known, whereas it is difficult to estimate 
this quantity when Ha emission line is not measured (which is 
the case when one wants to used [On] A 3727 instead to derive 
SFR). 

The paper is organized as follows: we first present the sam- 
ple we have used and the selection we applied to it (Sect. [2]), then 
we calibrate the SFR as a function of emission-line luminosities 
when a correction for dust attenuation is available (Sect. O, or 
when it is not available (Sect.|4]). We finally draw some conclu- 
sions (Sect. [5]). 



2. Description of the sample 

The physical properties of a total of 567486 galaxy spectra 
inside SDSS Data Release 4 (DR4, lAdelman-M cCarthv et al. 
l2QQ6l) have been made publicy available on the following web- 



site: jhttpjy/www^mpa-^garching . mpg.de/SDSS/[ Taking ad- 
vantage of this huge amount of data, we made a sub catalog with 
the SFR, the emission-line measurements, the spectral classifi- 
cation and the metallicities of a sub sample of galaxies. The se- 
lection will be described later in this section as we first describe, 
for the benefit of the reader, the main ingredients of the methods 
used to derive the physical parameters. 



2.1. Physical properties 

The spectral properties o f SDSS galaxie s have been measured 
with "platefit" software by ITremonti et al.l (|2004|) . For each spec- 
trum, they perform a careful subtraction of the stellar continuum 
and absorption lines by fitting a linear co mbination of single 
stella r population models of different ages (Bruzual & C harlotI 
l2003l) . Then they fit all emission lines at the same time as a 
unique nebular spectrum in which each line is modeled by a 
Gaussian. Thanks to the subtraction of the stellar component. 



B aimer emission lines are automatically corrected for underly- 
ing absorption. We note that the [Oil] A 3727 line that we use in 
this study is the sum of two Gaussians of the same width mod- 
eling the [On] A A 3726,3729 doublet. The AGN classification is 
based on the [Oiii]A5007/Hj3 vs. rNiilA6584/ Ha emission-line 
diagnostic diagram and has been performed by lKauffmann et all 
(2003ai) . 

The SFR and gas-phase oxygen abund ances have been com- 
puted from emission-lines fluxes usin^ the [Chariot & Longhettil 
(2001, hereafter CLOl) method bv lBrinchmann et al.1 (12004') and 
[Tremonti et al.1 (|2004|) respectively. This method compares all 
emission-lines fluxes together to a set of theoretical nebular 
spectra, which are modeled with five parameters: the metallicity, 
the ionization degree, the dust-to-metal ratio, the dust attenua- 
tion, and the SFR efficiency factor. For each galaxy, one estimate 
of each obs erved parameter is computed through a Bayesian ap- 
proach (see iBrinchmann etani2004l for more details). We note 
that the SFRs are computed assuming the iKroupal (120011) ini- 
tial mass function (IMF). One may s cale th e SFR estimates dis- 
cussed in this paper to the Salpete? ("1955") IMF by multiplying 
them by a factor 1.5 (or by adding 0. 176 dex to their logarithm), 
and to the lChabriej (l2003b IMF by multiplying them by a factor 
0.88 (or by subtracting 0.056 dex to their logarithm). The CLOl 
SFR will be used as the reference SFR throughout this study. 

We compute the luminosity L of each line starting from its 
measured flux F, and the redshift z of the galaxy, using the fol- 
lowing equation: 



L = 4;r(4-(l+z)-/o/(^0-^/W)'xF 



(1) 



The luminosities are calculated with the same cosmology 
ba sed on WMAP results (jSpergel et aP l2003l) as the one used 
by IBrinchmann et al.l (l2004 to estimate the SFR: Hq = 70 km 
s"^ Mpc~^ Qa = 0.7 and 0,^ = 0.3. We note that aU the caH- 
brations provided in this study are independent of the cosmology. 
They indeed compare SFR and emission-line luminosities which 
are affected in the same way by the distance modulus. 

Finally, w e have cross-m atched this catalog with the VAGC 
catalog (B lanton et al .120051) to get the ^-corrected absolute mag- 
nitudes which are used in this study. The absolute magnitude in 
the 5-band is derived from SDSS u^- and g'-bands using Eq. 1 of 
Moustakas et al. (20061 



2.2. Dust attenuation 

We will use the following notations in all this work: the ob- 
served flux of the three emission lines used in this study are de- 
signed by the name of the line: Ha, Hj3 and [Oil] (stands for 
[On] A 3727). The intrinsic flux corrected for dust attenuation is 
designed by a superscript /: Ha\ Hj3^ and [Oll]\ The attenuation 
law is designed by the function t(A), and the effective dust atten- 
uation in V-band by the symbol Ty. For simplification purposes, 
we also note: t{V) = t(5500A) (V-band), t(j3) = t(4861A), 
T(a) = t(6563A), and t(Oii) = t(3727A). The intrinsic flux 
f{{^) at any wavelength A as a function of the observed flux 
f^ (A) is given by the following equation: 



/i(A)=A(A)-exp 



T(y) 



.T(A) 



(2) 
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We compute the dust attenuation by comparing the observed 
to the intrinsic Ha/Hj3 emission-line ratios, using the following 
equation: 



.Baimer ^ ln(Ha/Hj3) - InfflaVHjgQ 



• T(y) 



(3) 



In the following work, we use the standard case B re- 
combination ratio (hereafter st andard intrinsic Baimer ratio): 
HaVH jS^ = 2.85 dOsterbrockl Il989h and the llCharlot & FalU 
(l200Qh mean attenuation curve: t(A) oc The choice of this 

curve is made to be consistent with the setup used to compute 
the SFRs in the SDSS DR4 data (the SFR is computed self- 
consistently with correction for dust attenuation). 

The dust attenuation recovered with an assumed constant in- 
trinsic Baimer ratio (ty^^^^^) does not take into account the vari- 
ations of this ratio with other physical parameters (e.g. metal- 
licity). We thus also use in this study the true dust attenu ation 
recovered with CLOl models by iBrinchmann et aP (|2004|) . We 
design it by the symbol Ty^^^ . The intrinsic flux corrected for 
dust attenuation using this parameter is designed by a superscript 
iC: Ha^c, Hj3^^ and [Oiif^. 

Finally, another way to estimate the dust attenuation is to 
compare the observed colors of the galaxies to a library of stellar 
population synthesis models (hereafter SED fitting). This work 
has been carried on with Bruzual & Chariot ( 2003) models on 
SDSS galaxies by Kauffmann et al. ( 2003b), who provide an 
estimate of the conti nuum attenuation in the z-band (A^) with 
ICharlot & Fall (l2000l) attenuation curve. It is straightforward to 
convert the parameter to Ty using the following formula: 



^SED 



A, /5500\ 
1.086' V8800 / 



-0.7 



1 



(4) 



Please note that the 1.086 factor is the conversion from mag- 
nitudes to opacities and A^/Ag is the mean ratio between the 
dust attenuation affecting stars (which dominate the colors of 
the galaxies) and the one affecting gas. The dust attenuation es- 
timated from SED fitting is characterized by a high uncertainty 
but might be used in a statistical way. 

2.3. Sample selection 

Thanks to the spe ctral classification provided by 
IKauffmann et al.l (l2003ah . we select only star-forming galaxies 
and remove all AGN from our sample. In the whole study, 
we additionally remove all galaxies for which Ha and Hj3 
emission lines (used for estimating the dust attenuation) are 
not measured with a signal-to-noise ratio of at least 5. We 
apply the same selection to the [On] emission line in order to 
derive the calibrations in which this line is involved. We note 
that the signal-to-noise ratios have been corrected according to 
the values provided in the SDSS/MPA website (based on the 
analysis of duplicated observations): they have been divided 
by 2.473 for Ha line, 1.882 for H^ line and 2.199 for [On] 
line. We also kept only galaxies with an available estimate of 



^CLOl 



. This selection was applied to remove a number of spectra 
with unreliable flux measurements. Please note that a number 
of objects may show a negative Ty^^^^^ estimate. This is due to 
statistical variations around the wrong assumption of a constant 
intrinsic Baimer ratio: these objects have not to be removed. 

Finally, some objects in SDSS DR4 catalog have been ob- 
served twice or more. We decided to keep only one observa- 
tion of each duplicated galaxy by selecting the one with the 
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Figure 1. This plot shows the histogram of the dust attenuation 
(left), or of the gas-phase oxygen abundance (right), both esti- 
mated with the CLOl method, for the star-forming galaxies in the 
SDSS DR4 data. The solid vertical line shows the mean value of 
the distribution, and the dashed vertical lines the l-O" tails. The 
values are given in the plot. The two dotted histograms shows 
the distribution of our two sub-samples (see text for details). 



best signal-to-noise ratio on the Ha emission line. We end up 
with 123 713 and 84733 star-forming galaxies for Ha-based and 
[On] -based studies respectively. 

All the parameters used in this study (SFR, metallicity, dust 
attenuation) are the median estimates. 



2.4. Properties of the studied sample 

Fig- [D (left) shows the distribution of the effective dust attenua- 
tion in the V-band (estimated with the CLOl method) in the sam- 
ple we have obtained. The mean value is 1 .21 and the dispersion 
is 0.637. We note that this values are obtained with the selection 
on Ha and Hj3 lines only. The new mean and dispersion when 
we add the selection on the [Oil] line are 1.10 and 0.547 respec- 
tively. Fig. [T] (right) shows the distribution of the gas-phase oxy- 
gen abundance in the same sample (for objects with an available 
measurement of this parameter). The mean value is 8.89 with a 
dispersion of 0.24 (8.88 mean value with a dispersion of 0.21 if 
the selection on the [Oil] line is applied). 

Fig. [T] also shows the distributions of four sub-samples (two 
in each panel). These sub- samples have been randomly selected 
from the main sample, forcing their mean values to be closest as 
possible to the mean value of the main sample plus or minus its 
dispersion, and their dispersion to be half the dispersion of the 
main sample. We end up with four sub- samples: low dust atten- 
uation, high dust attenuation, low metallicity, and high metallic- 
ity. They will be used to test the calibrations on samples with 
different dust properties. 



3. SFR calibration with an available dust estimate 

We present in this section emission-line calibrations of the SFR 
based on dat a alr eady corrected for dust attenuation. We first fo- 
cus in Sect. 13.21 on data with dust attenuation estimated from 
the Baimer decrement, which is the common case in the liter- 
ature. We then study in Sect. 13.31 the different calibrations we 
obtain when using the dust attenuation estimated from the CLOl 
method or SED fitting. 



4 



B. Argence and R Lamareille: Emission-lines calibrations of the Star Formation Rate from the Sloan Digital Sky Survey 



3.1. Base theory 

The common way to calibrate the SFR against one emission line 
luminosity (e.g. Ha) is a simple scaling relation of the following 
form: 



SFR = L(Ha)/r7Ha 



(5) 



This formula, where rfna is the efficiency factor, repro- 
duces the simple relation between the amount of ionizing pho- 
tons (the SFR) and the emission-line luminosity. However, it has 
been already shown that the Ha efficien cy factor depends on 
other physical properties of the galaxies: Chariot & Longhetti 
(l2QQll) have shown that r}Ha depends on the metallicity, metal- 
rich stars being less luminous and thus producing less Ha pho- 
tons in the interstellar m edium for a same SFR. More recently, 
iBrinchmann et al.l (|2QQ4|) have shown that rfna is also linked to 
the stellar mass of SDSS galaxies, massive galaxies producing 
less Ha photons than dwarfs for a given SFR. We note that these 
two results are obviously linked by the ma ss-metallicity rela- 
tion (iTremonti et al.ll2004l: iLamareille et al.h20Q81 Fortunately, 
IBrinchmann et al.l ( l2004l) have also shown that this effect is can- 
celed by the commonly used wrong assumption of a constant 
intrinsic Ha/Hj3 intrinsic B aimer ratio. In fact this ratio in- 
creases with metallicity, ending up with an overestimate of the 
dust attenuation for metal-rich galaxies. This latter effect almost 
exactly compensates the decrease of rjHa when metallicity in- 
creases. 

Nevertheless, we decided to adopt a more general approach 
using a power-law as described by the following formula: 



log(SFR) = eHalogL(Ha) -log(7]Ha) 



(6) 



The parameter is the exponent of the power law. The 
simple calibration described by Eq.[5]corresponds to the Eua = 1 
special case. We show in the following sections that this more 
general approach is not mandatory for the cases with a dust 
attenuation estimated from the B aimer decrement (^Ha is very 
close to 1), but that it is useful in other cases. 

We note that the correlation between dust attenu ation and 
metallicity ( Cortese et ani2006l: iMouhcine et al.ll2005l) also has 
an effect on the slope. However the goal of our work in not to 
study in details such phenomenon. 

In the following sections, the efficiency factors are expressed 
in erg s~^ (M© yr~^)~^ units and the exponents of the power 
laws are unit-less. All fits are least square fits with errors in x 
and y. The errors on the fitted parameters are given by the rms of 
40 bootstrap estimates. The term "dispersion" relates to the rms 
of the residuals around the fitted calibrations, the term "shift" to 
the mean of these residuals. The given residual slopes are unit- 
less. 



3.2. Dust estimated from the Balmer decrement 

We now present possible calibrations of the SFR based on data 
corrected with dust estimated from the Balmer decrement (see 
Eq. [3]). We remind the reader that the results presented in this 
subsection are valid only if we correct the observed flux for dust 
attenuation using Eqs.[2]and[3l i.e. making the wrong assumption 
that the intrinsic Balmer ratio is constant over all the range of 
stellar masses. 



3.2.1. Improved standard calibrations 

iKennicuti (|l998) has developed two calibrations of the SFR ver- 
sus Ha^ or the [Oll]^ emission lines which are now widely used 
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mean=0.102 
rms = 0.178 
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Figure 2. These four plots show the residuals of the SFR recov- 
ered from an emission-line calibration as compared to the ref- 
erence CLOl SFR, for star-forming galaxies in the SDSS DR4 
data. This comparison is shown as the logarithm of the ratio 
between the derived and the reference SFR, as a function of 
the logarithm of the reference SFR (in M0/yr). The number of 
objects used in the comparison, the shift, the dispersion, and 
the Spearman rank correlation coefficient are given in the plot. 
The dashed horizontal line is the y = curve, the solid line 
is the y = a ■ x-\-b curve where a and b are the parameters 
of the linear regression given in the plot. Isodensity contours 
are overplotted i n white. The st udied calibration s are: top-left: 
Kennicutt ('1998') Ha^; top-right: iKennicuttl (1 1 998h (corrected by 
Kewley et al. 2004) [Oii]^; bottom-left: our improved Ha^ cafi- 
bration; bottom-right: our improved [Oll]^ calibration. 



in the literature. We take advantage of the SDSS DR4 data to test 
these two calibrations as compared to the high quality CLOl esti- 
mate of the SFR. We also derive two new improved calibrations 
based on Eq. [6] rather than Eq.[5] 

Fig. [2] (top-left) shows the comparison b etween the CLOl 
SFR and the one recovered with the standard [kennicutti (Il998h 
Ha^ calibration, which uses log(7]y^) = 41.28 if we scale their 
results to our adopted IMF. We find a good agreement: the dis- 
persion is 0.17 dex, which is very similar to the minimum un- 
certainty associated to the CLOl SFR estimates (i.e. 0.16 dex). 
Thus, we confirm that the Ha emission Une luminosity alone is 
sufficient to recover almost perfectly the SFR. 

However, we see that the reduced efficiency factor and the 
dust attenuation overestimate do not exactly cancel each other 
on the massive end of the sample, which induces a mean of the 
residuals of +0.1 de x. One might decide to subtract 0.1 dex to 
the Kennicutti (119981) calibration in order to find a null mean of 
the residuals, or we can try to derive an improved calibration. 

Fig. [3] (left) shows our improved calibration for the correla- 
tion between the Ha^ emission-line luminosity and the SFR in 
SDSS DR4 data. We find the following best fit values: 



39.38 ±0.03 
0.951 ±0.001 



(7) 
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log L(HaO (erg/s) log L([OII]A37270 (erg/s) 

Figure 3. Relation between the SFR (logarithm of Mo/yr) and the Ha^ (left) or [Oll]^ (right) emission lines luminosities (logarithm 
of erg/s) for star-forming galaxies in the SDSS DR4 data. The solid line is a least-square fit to the d ata (errors on x and y). The 
dashed line is the Kennicutt (1998) calibration corrected to our assumed IMF, and corrected for dust by .Kewley et alJ (|2QQ4|) on the 
right panel. Isodensity contours are overplotted in white. 



One can notice that the formal errors of these two parame- 
ters are very low thanks to the large number of objects, which 
makes this calibration quite reliable. The uncertainty on the ex- 
ponent of the power law especially tells us that this parameter 
is very well constrained, while not equal to 1 as compared to all 
previous studies. As shown in Fig. [2] (bottom-left), the mean on 
the residuals is now almost null (0.01 dex) with the improved 
calibration. A small residual slope (—0.04) is present as the im- 
proved calibration is now less accurate for very low SFR values. 
This residual slope is not significant given the low Spearman 
rank correlation coefficient: —0.23. 

We now consider the case where the Ha line is not observed, 
or not measurable, and that we decide to use the [On] A 3727 
emission line instead. In such case, the dust attenuation might 
be derived from another B aimer ratio not involving the Ha line 
(e.g. H7/Hj3,H5/H7, ...) 

As already stated in the literature, the relation between the 
SFR and the [Oil] line has to be considered with caution because 
of the stronger dependence of this line on metallicity (com- 
pared to Ha), as well as a stronger sensibility to dust attenua- 
tion. However, as shown in Fig. [2l (top-righ t), we find a good 
agreement with the standard Kenni cutt (11998 ) [Oll]^ calibration, 
which uses log(T]^Qjj^ )= 41 .36 if we correct their results for dust 
attenuation (iKewley et al .112004)) and scale them to our adopted 
IMF. 

It seems that the effect of metallicity on the [Oil] line does 
not produce any systematic shift and only leads to a higher dis- 
persion: 0.25 dex. To be more precise, the three effects of the 
metallicity on the relative strength of the [Oil] line, on the effi- 
ciency factor, and on the estimate of the dust attenuation exactly 
compensate in order to produce a exponent of the power law very 
close to unity. 

Fig. [3] (right) shows our improved calibration in the correla- 
tion between the [Oll]^ emission-line luminosity and the SFR in 
SDSS DR4 data. We find the following best-fit values: 



logT]^ = 41.94 ±0.7 



[On] 



1.014±0.02 



(8) 
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Figure 4. Relation between the [Oii]A3727VHa^ emission line 
ratio and gas-phase oxygen abundance (estimated with the CLOl 
method) for the star-forming galaxies in the SDSS DR4 catalog. 
The dashed curve is a semi-empirical estimation of this relation 
(iKewlev et al.ll2004l) . 



As expected, and as shown in Fig. [2] (bottom-right), this cali- 
bration does not produce any significant improvement compared 
to the iKennicuttI (11998 ) (corrected by Kewley et al.l l2004D cal- 
ibration. In both cases, the observed residual slope of —0.09 is 
not significant (Spearman rank correlation coefficient of —0.24). 

3.2.2. New calibrations 



As already s hown in previous works ("Kewley et alj 120041: 
iMoustakas et al. 2006; Weiner et al. 2007), the best way to pro- 
vide a reliable SFR vs. [On] calibration is to actually use the 
SFR vs. Ha^ calibration, and to calibrate either the [Oll]/Ha^ 
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Fig ure 5. Same le gend as in Fig.O The studied calibrations are: 
left: lKewley et al.l(l2QQ4 [Oii]^ with a correction for metalHcity; 
right: our new [Oll]^+Hj8^ calibration. 

(observed-to-intrinsic) or the [Oll]VHa^ (intrinsic-to-intrinsic) 
line ratio against a well chosen parameter. Then the star forma- 
tion rate is given by the following formula: 

logSFR = e^i^ (logL([Oii]^-)-log([Oii]iVHa^)) -\ognk^(9) 

iKewlev et aP (l2QQ4l) has calibrated the [Oii]VHa^ line ra- 
tio against the gas-phase oxygen abundance (12 + log(0/H)), 
ending up in a n appreciable in i pro vement of the calibration. As 
also shown by Mouh cine et al. ^ (2005) with 2dFGRS data and in 
Fig.m there is indeed a robust correlation between these two pa- 
rameters. The dashe d line in Fig [ 1] show s the relation derived 
semi-empirically by ' Kewley et al ] (12004 . The residuals around 
this relation are characterized by a mean of —0.04 dex and a rms 
of 0. 1 dex in SDSS DR4 data. 

However, we see three drawbacks against the use of this re- 
lation: 

- First the gas-phase oxygen abundance is not an easy pa- 
rameter to derive. There are many different calibrations 
whic h can be different from each other by up to 0.5 dex 
(see iKewlev & EllisonI l2008l for a detailed discussion). 
Moreover, most of these calibrations need many emission 
lines to be measured and we emphasize that, in such cases, 
it might be easier to use directly the SFR estimated from 
multiple-lines fitting methods like CLOl. 

- Second, we already stated before that correcting the [Oil] 
emission line for dust attenuation requires in most common 
cases that the Ha emission line is being observed: in this 
case one would estimate the SFR directly from the Ha line, 
rather from the [Oil] line. Thus, trying to improve the [Oil] 
calibration with data corrected for dust attenuation is actually 
not applicable in most common cases. 

- For the other cases (e.g. dust attenuation estimated from an- 
other B aimer ratio), we also remark that using the metallicity 
estimate (or any other derived parameter) to correct the cali- 
bration is not necessary. Indeed a very simple relation, which 
do not need a metallicity estimate, already exists between 
[Oll]\ Ha^ and Hj3^ emission lines. This simple relation is 
expressed by the following equation assuming the standard 
intrinsic B aimer ratio: 



M 

Hai 



[On]i tW _ [Onf 



(10) 



which an estimate of the gas-phase oxygen abundance was avail- 
able. Before doing that, we have corrected our available CLOl 
metallicities to Kewlev & Dopita ( 2002) metallicities, using the 
formula in Table B3 of iKewlev & ElHson (20 08|). 

As shown in Fig. [5] (left), the Kewley et al.i (12004 ^ calibration 
shows a non-negligible 0.1 dex shift. Moreover, the dispersion 
of the relation is 0.26 dex, which is not better than the standard 
[Oll]^ calibration without a correction for metallicity. Their is 
also a small, not significant, residual slope of —0.09 (Spearman 
rank correlation coefficient of —0.23) . 

Eqs. [9]and[T0lare now combined to provide an estimate of 
the SFR, which is compared to the CLOl SFR on Fig. [5] (right). 
The relation is, as expected, very good as no dispersion is added 
by Eq.[TOl The small dispersion of the relation (0.16 dex) comes 
only from the uncertainty of the Ha^ calibration that we derived 
previously. This calibration do not show a significant residual 
slope (Spearman rank correlation coefficient: —0.17). In fact this 
new [Oii]^+Hj3^ calibration is expected to give exactly the same 
results as the underlying Ha^ calibration since the relation be- 
tween them is the constant 2.85 factor. Fig. [5] (right) is nothing 
else than Fig. [2l (bottom-left) without some missing points be- 
cause of the additional selection on the [Oil] line. These missing 
points explain the smaller scatter. 

Following this conclusion, we emphasize that a simple H/3^ 
calibration is even easier to derive and to apply to the observa- 
tions using this relation: 



Ha^ = Hj3^x2.85 



(11) 



We applied Eq. 15 of i Kewley et al.l (|2004|) in order to compare 
their metallicity- and dust-corrected calibration to the reference 
CLOl SFR. We have done that only for the 84730 galaxies for 



The Hj3Ms expected to give exactly the same results as the 
Ha^ calibration, provided that a correction for dust attenuation 
has been estimated from the assumption of a constant intrinsic 
B aimer ratio. 

3.3. Dust estimated from another method 

The dust attenuation might not be estimated with the assump- 
tion of a constant intrinsic B aimer decrement. Other methods 
such as CLOl can provide a better estimate of the dust attenua- 
tion, not biased towards the metallicity dependence of the intrin- 
sic B aimer ratio (thanks to the use of all available emission-line 
measurements). The dust attenuation estimated from SED fitting 
methods does not either make the assumption of a constant in- 
trinsic B aimer ratio. 

As stated before, the calibrations derived in previous subsec- 
tion should show an exponent of the power law greater than 1 .0, 
because of the metallicity dependence of the relation between 
the SFR and the emission-line luminosities (assuming that the 
higher is the SFR the higher are the stellar mass and the metal- 
licity). However, this effect is not seen since it is almost exactly 
compensated by the overestimate of the dust attenuation at high 
metallicities. 

We thus study in this subsection the results obtained with 
dust estimated from other methods. 

3.3.1. Metallicity-unbiased calibrations 

We provide new calibrations to be used on data corrected with 
a good estimate of the dust attenuation, i.e. not biased towards 
metallicity. To this end, we corrected the emission-line luminosi- 
ties using the dust attenuation estimated with the CLOl method, 
instead of the one estimated from the B aimer decrement. We em- 
phasize that this work is mainly done for comparison purposes: 
using the CLOl method to estimate the dust attenuation makes 
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Figure 6. Relation between the SFR (logarithm of Mo/yr) and the Ha^^ (left), [Oll]^^ (center) or Hj3^^ emission lines luminosities 
(logarithm of erg/s) for star-forming galaxies in the SDSS DR4 data. The solid line is a least-square fit to the data (errors on x 
and y). These data have been corrected using the unbiased CLOl estimate of dust attenuation. The dashed lines show the previous 
calibrations obtained with dust estimated from the assumption of a constant intrinsic Balmer ratio (see Fig.O. Isodensity contours 
are overplotted in white. 



our SFR calibration useless since the CLOl SFR is already better 
constrained. 

Fig. [6] shows the relations between the SFR and Ha^^, 
[Oll]^^ or Hj3^^ emission lines. In the three cases, we find as 
expected an exponent of the power law greater than 1 . However 
we see that the metallicity dependence of the efficiency factor 
has not a strong effect on the slope for the Ha^^ and Hj3^^ cali- 
brations. We remind the reader that, in contrary to what we have 
said in Sect. 13.21 the calibrations based on Ha^^ or Hj3^^ lines 
are now slightly different since the ratio between these two lines 
is not any more assumed constant. We obtain the following re- 
sults: 



^iC 



iC 
a 



logT] 



r.iC 

-[Oil] 



iC 
[On] 



log 7]^ 

ciC 



iC 

m 



44.22 ±0.03 
1.072 ±0.001 

= 49.90 ±0.9 
= 1.213±0.03 

44.60 ±0.03 
1.094 ±0.001 



(12) 



(13) 



(14) 



In the three cases, the dispersion is reduced thanks to the 
better estimate of the dust attenuation provided by the CLOl 
method. Is it now < 0.1 dex in the Ha^^ and Hj3^^ calibrations 
and 0.24 dex in the [Oll]^^ calibration. 



3.3.2. Dust estimated from SED fitting 

One may ask if emission lines can also be used to derived SFR 
with an estimation of dust attenuation coming from SED fit- 
ting. Fig. [7] (left) shows the correlation between the intrinsic 
Ha emission-line luminosity and the SFR in SDSS DR4 data, if 
we correct Ha lumin osities with the dust at tenuation computed 
from SED fitting by KauffinaimetaD ([2003b) (see Eq.ll. Note 
that only the 87213 galaxies with an available estimate of the 
dust attenuation are plotted. 

We use a stellar- to-gas attenuation ratio of 0.3 which 
corresponds to the mean value observed in the SDSS data 
(Brinchmann, private communication). We see that the exponent 
of the power law is quite smaller (0.710 ±0.001), telling us that 
the SED fitting dust attenuation is even more overestimated at 
higher masses than with the Balmer decrement method. 




log L(Ha's™) (erg/s) 



log L(Ha's™) (erg/s) 



Figure 7. Relation between the SFR (logarithm of Mo/yr) and 
the Ha^^^^ emission lines luminosity (logarithm of erg/s) for 
star-forming galaxies in the SDSS DR4 data. The solid line is 
a least- square fit to the data (errors on x and y). The data has 
been corrected for dust attenuation estimated with the SED fit- 
ting method (see Eq.|4l). The stellar- to-gas attenuation ratio has 
been assumed to 0.3 (left: SDSS mean value), or 0.44 (right: 
FCalzetti 2001). The dashed line is the Kennicutt (1998) calibra- 
tion corrected to our assumed IMF. Isodensity contours are over- 
plotted in white. 



But no strong conclusion can be drawn: the dispersion that 
we find (0.40 dex) is indeed quite high and moreover, it does 
not take into account the likely variations of the stellar-to-gas 
attenuation ratio between the SDSS data and any other sample. 
Using instead a stellar-to-gas attenuation ratio of 0.44 (Calzett^ 
[2001), we found in Fig. [7] (right) an exponent of the power law 
of 0.864 ±0.001. The derived efficiency factors with a mean 
stellar-to-gas attenuation ratio of 0.3 and 0.44 are respectively 
29.67 ± 0.02 and 35.79 ± 0.03. We emphasize that a stellar-to- 
gas attenuation ratio of 0.44 is not the right value representing 
these data, even if it gives by chance a less dispersed calibration 
(0.34 dex). 

We conclude that the most critical issue in using dust es- 
timated from SED fitting, in order to correct emission lines, 
comes from the uncertainty in the stellar-to-gas attenuation ra- 
tio. This problem makes SFR derived using this method almost 
completely random on a galaxy per galaxy basis. We note that 
this method might be used in a statistical way, but with great un- 
certainties coming both from the dispersion observed in Fig. [71 
and from the fact that one has to know the right stellar-to-gas 
attenuation ratio to apply to his data. 
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4. SFR calibration without a dust estimate 

We present in this section emission-line calibrations of the SFR 
based on data not corrected for dust attenuation. Because of the 
small wavelength coverage of many spectroscopic surveys, it is 
common that not all emission lines necessary to derive a correc- 
tion for dust attenuation, and to compute a SFR, are present in a 
given spectrum. 

There are two ways to handle such data: ( i) use an assumed 
mean extinction and apply a statistical correction for dust attenu- 
ation (see Sect. l4.1l) : or ( ii) use a SFR calibration which provides, 
still in a statistical approach, a self-consistent correction for dust- 
attenuation, i.e. with dust properties recovered from another pa- 
rameter (see Sect. 14.21) . We remark that such self-consistent cor- 
rection depends on the properties of the sample used to do the 
calibration. Thus, we present in Sect. 14.31 a way to correct for 
this bias. 



4.1. Use of an assumed mean correction 

When dust attenuation cannot be reliably derived from the data, 
it is common in the literature to use an assumed mean correction 
(frequently Ay = 1). 



4.1.1. Standard calibrations with Ay = I 

In this subsection, we discuss the quality of the standard calibra- 
tions used on SDSS DR4 with an assumed mean correction for 
dust of Ay = 1 . 

One important question we need to answer is: what kind of 
calibration should we use on data corrected with a mean dust 
attenuation ? Should we use the standard calibrations derived 
with the wrong assumption of a constant intrinsic B aimer ra- 
tio (Sect. 13.2b , or our new calibrations derived with the better 
CLOl estimates of the dust (Sect. l3.3l) ? The answer basically de- 
pends on the way the mean dust attenuation have been estimated 
(is it biased towards the metallicity dependence of the intrinsic 
B aimer ratio or not ?). 

We thus plot the comparison of both calibrations in Fig. [S] 
The main conclusion we draw from this figure is that the resid- 
uals are not only highly dispersed (rms in the range 0.3-0.5 dex 
depending on the calibration) but also show a clear residual slope 
in all cases (—0.2 to —0.5 depending on the calibration). This 
residual slope is the signature of the correlation between the 
dust attenuation and the SFR in galaxies. This correlation is ob- 
viously not taken into account by the assumption of a constant 
dust attenuation for the whole sample. 

We note that the residual slopes are significant as shown by 
the Spearman rank correlation coefficients (—0.5 to —0.7). It is 
clear also that the residual slopes obtained with the Ha^^, [Oll]^^ 
, and Hj3^^ calibrations, based on the CLOl estimates of the dust, 
are not only smaller but also less significant (smaller Spearman 
rank correlation coefficient). 

The method of the assumed mean correction thus leads to 
poor results, either on galaxy-per-galaxy studies (high disper- 
sion) or on statistical studies (non-negligible residual slope). 
This method applied with the [ Oil ] line is particularly inadvis- 
able (worst dispersion and steeper residual slope). 
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Figures. Same legend as in Fig. [2l The studied calibrations 
are the standard Ha^ (top-left), [Oii]^ (middle-left), and Hj3^ 
(bottom-left) calibrations, or the metallicity-unbiased Ha^^ (top- 
right), [Oll]^^ (middle-right), and Hj3^^ (bottom-right) calibra- 
tions, all applied on data corrected using an assumed mean dust 
attenuation Ay = 1. 



4.1.2. Additional sources of biases 

We now discuss the additional effect of the use of a wrong mean 
correction. Following Eq. [2] and [6l the bias introduced on the 
SFR can be computed like that: 



log 



/SFR^ 



V SFR 



£ X 



Ty 



In(lO) 



(15) 



where Ty is the wrong assumed dust attenuation, and SFR^ is 
the wrong derived SFR. Below, we investigate three possible 
sources of bias. 

The first bias comes from the assumption of Ay = 1 . The 
sample used in this study actually shows a mean dust attenuation 
^^CLOi^ = 1.21 (see Sect.[2]and Fig.[T]). Converting from opac- 
ities to magnitudes, we obtain (Ay^^^) = (t^^^^) x 1.086 = 
1.31. The difference between the assumed Ay = 1 and the actual 
value in the sample is partly responsible of the non-null means 
of the residuals observed in Fig. [H applying Eq. [151 we find a 
theoretical bias of the order of —0.15 dex, which is similar to 
the observed values for the CLOl calibrations (Fig. [8] right). 
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Figure 9. Relation between the dust attenuation estimated with 
the CLOl method or with the B aimer decrement method, for star- 
forming galaxies in the SDSS DR4 data. The dashed line is the 
y = X curve. The solid line is the y = a-x-\-b curve where a and 
b are the parameters of the least- square fit (errors in x and y), as 
given in the plot. Isodensity contours are overplotted in white. 



The second bias comes from the dispersion in the distribu- 
tion of dust attenuation. In the SDSS DR4 data, the dispersion on 
the derived dust attenuation is 0.64, which means that the mean 
dust attenuation may not be assumed with a better precision than 
±0.64. Thus it implies a possible systematic shift of the order of 
±0.30 dex. As stated in Sect. [2l the mean dust attenuation is 
^^CLOi ^ = 110 instead of 1 .21 when we add the selection on the 
[On] line: this implies a bias of the order of ±0.05 dex. 

Finally, we remind the reader that he has to think carefully of 
the origin of the assumed mean correction before applying it, i.e. 
he has to check whether this assumption comes from a dust at- 
tenuation derived with a constant B aimer decrement, or from an 
unbiased estimation. Having checked that, he would be able to 
apply the right calibrations. Fig. [9] gives for convenience the re- 
lation between t^^^"^^^ and the unbiased T^^^^ in the SDSS DR4 



data. We fit an empirical relation which follows this formula: 



Balmer^^^53^ CLOl. 



0.3 



(16) 



This implies a theoretical bias of the order of —0.25 dex, which 
is again similar to the observed values for the standard calibra- 
tions (Fig.[8]left). 

Hence we conclude that the use of an assumed mean correc- 
tion has to be handled with great care. 

4.2. Use of a self-consistent correction 

As stated in previous subsection, the use of an assumed mean 
correction has one main drawback: it does not take into account 
the variation of the dust attenuation as a function of SFR. We 
thus discuss in this subsection three possible ways to reduce this 
bias. The idea common to these three possibilities is that the 
relation between the dust attenuation and the SFR might be itself 
recovered from the data, i.e. self-consistently. 

We explore the self-consistent corrections obtained either 
from the observed line luminosities, another parameter like the 
magnitude, or a combination of two emission lines. 
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Figure 10. Same legend as in Fig. [2l The studied calibrations 
are: top-left: Kennicutt (1998) Ha^ applied on uncorrected data ; 
top-right: Kennicutt (1998) (corrected by iKewlev et aP l2004 
[Oii]^ on uncorrected data; middle-left: our new Ha calibration; 
middle-right: our new [Oil] calibration; bottom: our new Hj3 cal- 
ibration. 



4.2.1. Single-line calibrations 

The simplest way to apply a self-consistent correction for dust is 
to relate the dust attenuation to the line luminosity itself. Hence, 
it consists on calibrating the SFR directly with the observed, not 
intrinsic, line luminosity with a two-parameter law as described 
in Eq.[6l 

Fig. [T0| (top-left) sho^y s the results obtained when using the 
standard Kennicutt] (Il992l) Ha^ calibration on data uncorrected 
for dust attenuation. We clearly see that the dust attenuation 
causes an underestimate of the SFR (the mean of the residu- 
als is —0.48 dex), and that this effect is increasing with SFR. 
Moreover the dispersion is 0.32 dex, and there is a residual slope 
of —0.24 which is significant (Spearman rank correlation coeffi- 
cient of —0.63). Taking this into account, a two parameter law is 
necessary, differently from what we observed on dust-corrected 
data in Sect.[T2l 
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Figure 11. Relation between the SFR (logarithm of Mo/yr) and the Ha (left), [Oil] (center) or Hj3 (right) observed emission lines 
luminosities (logarithm of erg/s) for star-forming galaxies in the SDSS DR4 data. The solid line is a least-square fit to the data 
(errors on x and y). The dashed lines show the intrinsic calibrations obtained with dust estimated with the B aimer decrement method 
(see Fig. [3]). Isodensity contours are overplotted in white. 



Fig. [TT] (left) shows the correlation between the observed Ha 
emission-line luminosity and the SFR in SDSS DR4 data. We 
obtain the following new calibration: 



logT]Ha 



46.49 ±0.04 
1.141 ±0.001 



(17) 



As shown in Fig.[TOl(middle-left), the dispersion of this new cali- 
bration is 0.30 dex, which tells us that the SFR recovered without 
dust attenuation correction are almost twice more dispersed than 
when an estimate of dust attenuation is available. The mean of 
the residuals are now almost null (0.02 dex). This calibration is 
better, in terms of the dispersion, than the results obtained with 
an assumed mean correction (see Sect. 14.11 dispersion of 0.33 
dex) but, it still quite poor. 

Unfortunately the use of the two-parameters law, which al- 
lows in principle to take into account the dependence between 
the dust attenuation and the SFR, does not cancel completely the 
residual slope which is still of —0.13. However it is clear, from 
the comparison of Fig.[8]and Fig.[TOlthat this remaining residual 
slope is much less significant (Spearman rank correlation coeffi- 
cient of -0.36 instead of -0.5 to -0.7). 

No estimation of the dust attenuation would be available in 
most cases where one has to use the [On] or the Hj3 emission 
line in order to compute a SFR. This is basically due to the fact 
the Ha line is commonly used to estimate the dust attenuation, 
and that there is no reason to use another line if the Ha line is 
measured. 

Fig. [TOl (top-right) sho ws the results obtained when using the 
standard lkennicuttl (Il992h [Oll]^ calibration on data uncorrected 
for dust attenuation, the dispersion is 0.48 dex which is high. 
There is also a significant systematic shift of —0.83 dex, plus 
a significant (Spearman rank correlation coefficient of —0.68) 
residual slope of —0.47. Fig. [TT] (center) shows the relation be- 
tween the observed [On] emission line luminosity and the SFR 
in SDSS DR4 data. We obtain the following best-fit values: 



log^[Oii] 

^[Oii] 



60.06 ±1.9 
1.484 ±0.05 



(18) 



But, as we can see from Fig. [TO] (middle-right) the dispersion 
(0.54 dex) is higher than the results obtained with an assumed 
mean correction (see Sect. l4.2l dispersion of 0.48 dex): the corre- 
lation between observed [On] luminosity and SFR is rather poor. 



Nevertheless, this new calibration has the advantage of show- 
ing a smaller residual slope of —0.22, also much less significant 
(Spearman rank correlation coefficient of —0.25). 

Fig. [TT] (right) shows the relation between the observed Hj3 
emission line luminosity and the SFR in SDSS DR4 data, which 
gives the following parameters: 



logT7H/3 
%/3 



46.54 ±0.04 
1.159±0.001 



(19) 



As shown in Fig. [TO] (bottom), this new calibration shows a 
small improvement as compared to the method of an assumed 
mean correction: dispersion of 0.38 dex instead of 0.39 dex, 
smaller residual slope of —0.19, smaller Spearman rank corre- 
lation coefficient of —0.41 for this residual slope. 

The two-parameters calibrations presented in this subsec- 
tion do not show significant residual slopes as compared to the 
method of the assumed mean correction. Unfortunately, the gain 
of having no more significant residual slopes is diminished by 
the still high dispersion. We claim that these new self-consistent 
single-line calibrations are however useful when one as no idea 
of the mean dust attenuation he wants to assume. As stated be- 
fore, the use of wrong mean correction could lead to a non- 
negligible systematic bias. 

4.2.2. Discussion of previous works 

Another way to correct self-consistently for dust attenuation, as 
long as for other parameters (e.g. metallicity, ionization degree, 
...), is to use the absolute magnitude. Indeed, the absolute mag- 
nitude provides an estimate of the dust attenuation through the 
general correlation between mass and dust attenuation in galax- 
ies. It also provides an estimate of the metallicity through the 
luminosity-metallicity relation. We note that these calibrations 
are only derived for the [Oil] line. Again, it is better to correct 
directly for dust attenuation using the B aimer decrement when 

Ha is o bserved. 

The iMoustakas et al.l (l2006l) [On] calibration is based on 
a correction with 5-band /^-corrected absolute magnitude, and 
starts now to be widely used in the literature. However, since the 
correction is provided by the authors in a few number of discrete 
points (see Table 2 of their paper), it is not clear how to inter- 
polate between them. The common approach now used in other 
works is to calculate a global linear regression. Fig. [TJ] shows 
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Figure 12. Relation between the SFR efficiency factor of the 
[On] emission line and the 5-band ^-corrected absolute mag- 
nitude (Vega sy stem). The data po ints are the values given in 
Table 2 of Moustakas et alJ (l2QQ6h . The solid curve is our 2nd 
degree polynomial least square fit w^ith the coefficients given in 
the text. 
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Figure 1 3. Same legend a s in Fi g. O The studied calibrations 
are: left: [ Moustakas et al.1 (l2006h [OiiJ+M^ calibration; right: 
IWeiner et al.l (l2007h rOlll+MH^ calibration. 



The IWeiner et"aD (l2QQ7l) [On] calibration is based on 
a correction with ^-band ^-corrected absolute magnitude. 
The authors give the slope of the relation between the 
log([Oll]^^^/Ha^) ratio and the if -band ^-corrected absolute 
magnitude: +0.23 dex/mag. The zero-point is obtained by ap- 



that this approach is not valid in the whole domain. We thus pro- 
vide a parabolic fit which leads to the following final formula: 

log SFR = logL( [Oil] ) - 4 1 + 6. 86 + 0. 862 X Mfi + 0.027 X m|(20) 

where Mb is the 5-band ^-corrected absolute magnitude in the 
Vega system (^vega = ^ab + 0.09). 

Fig. [Is] (left) shows how the results obtained with the 

iMoustakas et al.i (2006) calibration compares with the reference 
CLOl SFR. This calibration shows a significant systematic shift 
of —0.22 dex. However, the dispersion of 0.43 dex is better than 
previous calibrations on the observed [Oil] luminosity. We ob- 
serve also a non-negligible residual slope of —0.30 (Spearman 
rank correlation coefficient of —0.46), but not stronger than the 
residual slope obtained with the method of the assumed mean 
correction. 



plying the 0.68 dex A(Ha) extinction to the log([Oll]^^^/Hao^^) 
ratio of —0.32 dex at Mh = —21. The final formula is: 

log SFR = 4^ (logL([Oii]) - 0.23 x - 3.83) - log t]^!^ (21) 

where Mh is the //-band ^-corrected absolute magnitude in the 
AB system. 

Fig. fT3l (r ight) shows how the results obtained with the 
IWeiner et aP (l2007 l) calibration compares with the reference 
CLOl SFR. The mean of the residuals is —0.24 dex, the dis- 
persion of 0.41 dex is better than the calibration on the ob- 
served [On] luminosity, but is still significantly high. Moreover, 
this smaller dispersion may be only due to the small number of 
galaxies (13678) with an available measurement of Mh in our 
sample. Finally, their is also a significant residual slope in this 
calibration of —0.36 (Spearman rank correlation coefficient of 
-0.55). 

We conclu d e tha t the calibration s prop osed by 
IMoustakas et al.1 (l2006 h and IWeiner et af] (120071) lead to 
poor results on the SDSS DR4 data. They do neither signifi- 
cantly better, nor significantly worse, than the results obtained 
with the method of the assumed mean correction, or with 
our direct calibrations between the SFR and observed line 
luminosities. 



4.2.3. Two-lines calibrations 

Finally we consider the general ca se where no est i mation 
of the d ust atten u ation is available. "Moustakas et alj (l2006h 
and Wei ner et"an (l2007l) have produced calibrations of the 
[Oll]/Ha^ line ratio against respectively 5-band or //-band ab- 
solute magnitudes. Following our idea of using only emission- 
line measurements, we explore a possible correlation between 
[Oll]/Ha^ (observed-to-intrinsic) and [Oll]/Hj3 (observed-to- 
observed) line ratios. Fig. [141 (left) shows the best-fit relation, 
as defined by the following equation: 

log([Oii]/Hai) = 3.243(±0.3)-log([Oii]/H/3)-1.74(±0.09) (22) 

The mean and the rms of the residuals around the fitted line 
are respectively —6 • 10~^^ dex and 0.37 dex. We note that the 
objects plotted upward the dashed line in Fig. [14] (left) have neg- 
ative dust attenuation when we derive them from Ha/Hj3 ratio 
with the assumption of a constant intrinsic B aimer ratio. 

Putting Eqs. [9] and [22] together, we obtain the following cali- 
bration formula: 

logSFR= (logL([Oii]) - 3.243 • log ([Oii]/H/3) + 1.74) - log Tlka(^3) 

Fig. [141 (center) shows the quality of the recovered SFR as 
compared to the reference CLOl SFR. We see that our new cal- 
ibration show no significant systematic shift (0.02 dex), and no 
significant residual slope (0.01 with a Spearman rank correla- 
tion coefficient of 0.06). Moreover our new calibration shows a 
dispersion of 0.29 dex which, while being still higher than the 
results obtained with a correction for dust attenuation, is the best 
result obtained in Sect.|4]between all other possibilities: assumed 
mean dust correction, direct calibration with observed line lumi- 
nosities, or self-consistent correction with the absolute magni- 
tude. 

We also explore the results obtained by doing a direct least- 
square fit of the SFR as a function of the [Oil] line luminos- 
ity and the [Oll]/Hj3 line ratio, without assuming the underlying 
Ha^ calibration. We find the following second-degree best-fit re- 
lation: 



log SFR = -39.21+0.982 • logL([Oii]) - 2.396 • log ([Oii]/Hj8) 

-0.3354- log ([Oii]/H/3)^ 



(24) 
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Figure 14. Left: Relation between the [Oll]/Ha^ (observed-to-intrinsic) and the [Oll]/Hj3 (observed-to-observed) line ratios for 
star-forming galaxies in the SDSS DR4 data. The solid line is a least-square fit to the data (errors on x and y). The dashed line 
is the Ha = 2.85 • Hj3 line. Isodensity contours are overplotted in white. Center and right: Same legend as in Fig.[2l The studied 
calibration are the new [Oll]/Hj3 calibrations (center: Eq.[23l right: Eq. ). 



Table 1. Coefficients to be used to correct the SFR derived from 
an observed calibration, as a function of the dust attenuation fol- 
lowing Eq. [25l The coefficients are given for a dust attenuation 
estimated with the CLOl or the Balmer decrement method. 



method 


CLOl 


Balmer 


calibration 


a 


b 


a 


b 


Hoc 


-0.45 


0.56 


-0.28 


0.44 


[On] 


-0.80 


0.88 


-0.51 


0.75 


Hi8 


-0.57 


0.70 


-0.37 


0.58 


[Oll]/Hj8 (Eq.ES 


-0.18 


0.22 


-0.082 


0.14 


[Oll]/Hj3 (Eq.IJl 


-0.34 


0.37 


-0.20 


0.30 



Fig. [14] (right) shows the quality of this other calibration as 
compared to the reference CLOl SFR. As compared to Eq. [23l 
the SFR derived with Eq. [24l show a null systematic shift and a 
smaller dispersion of 0.23 dex. Nevertheless this smaller disper- 
sion is diminished by a residual slope of —0. 12 (Spearman rank 
correlation coefficient of —0.31). 

Because it is much easier to use and gives better results, our 
calibrations based on only two observed line fluxes ([Oil] and 
Hj3) should be preferred in most cases with no available correc- 
tion for dust attenuation. The minimum uncertainty when dealing 
with data not corrected for dust attenuation is then only 0.23 dex. 
We discuss in the following section the different cases where it 
is better to use Eq.[23]or Eq.[24l 



4.3. The self-consistent correction on samples with different 
properties 

It is clear at this stage of our study that all th e calibrations de- 
rive d or studied in Sect . [4] (inc luding those by iMoustakas et al.l 
|2006 and Weiner et al.l l2007h . which are based on observed 
quantities, depend on the properties of the reference sample. 
Nevertheless, we emphasize that such statistically significant 
sample such as the SDSS DR4 data, with a large selection func- 
tion, is not expected to be strongly biased in favor of one partic- 
ular population of galaxies. 

It might be useful however to know how these calibrations 
are sensible to the variations of the properties of the studied sam- 
ple, especially in terms of dust or metallicity. 



a=0.080 b = 0.130 
spear = 0.223 
mean=0.089 
-rms = 0.262 




"20611 star-forming 
S/N S 5 

a=0.090 b = -0.111 
spear=0.182 
mean= -0.086 
-rms = 0.309 




log SFR (Mg/yr) 



a=-0.014 b=0.140 
spear=-0.039 
mean=0.147 
-rms = 0. 156 



log SFR (Mg/yr) 



a=-0.020 b = -0.214 
spear=-0.043 
mean=-0.219 
-rms = 0.210 




log SFR (Mg/yr) 



log SFR (Mg/yr) 



Figure 16. Same legend as in Fig. O The studied calibration is 
the new [Oll]/Hj3 calibration given by Eq. [23] (top) or Eq. [24l 
(bottom), applied on the two dust sub-samples defined in Sect.[2l 
Left: (t^^^^) = 0.63, right: (t^^^^) = 1.84. 



4.3.1. Dependence on dust 

We now study in Fig.[T5]the residuals of our new calibrations, as 
compared to the reference CLOl SFR, as a function of the dust 
attenuation estimated either with the CLOl method, or with the 
assumption of a constant intrinsic Balmer ratio. This has been 
done for the poor quality single-line calibrations (Ha, [Oil], and 
Hj8, see Fig.[T0|, and for the better quality [Oil]/ Hj3 calibration 
(see Fig. [14] center and right). 

In all cases, the means of the residuals are very low, which is 
expected as we always derived the best-fit calibration. 

For the single-line calibrations, we obtain very significant 
correlations between the residuals and the dust attenuation, with 
Spearman rank correlation coefficients of the order of —0.75 to 
—0.95. Thus, we can derive linear relations which may be used 
to correct these calibrations for a desired studied sample, for 
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Figure 15. Same legend as in Fig. [21 except that the residuals of the SFR calibration are now plotted as a function of the dust 
attenuation estimated with the CLOl method (top) or with the Balmer decrement method (bottom). The studied calibrations are, 
from left to right: Ha, [Oil], Hj3, [Oll]/Hj8 (Eq.[211), and [Oll]/Hj3 (Eq.[24|. All calibrations are related to observed quantities. 



which one has a rough estimate of the mean dust attenuation. 
The corrective formula as derived from Fig. [15] is: 



log (SFR"^"") = log (SFR) + a X Ty + 



(25) 



The a and b coefficients are summarized in Table [T] We note 
that this correction formula give a null correction for a dust at- 
tenuation equal to the mean value in the SDSS DR4 sample: 



^CLOl 



1.21 in the Ha and Hj3 cases, and for t\ 



€L01 



1.1 in the 



[On] case. 

For the [Oll]/Hj3 two-lines calibration given in Eq. [23l we 
see that the residuals do not show any more any strong signif- 
icant correlation with the dust attenuation: smaller slopes, and 
smaller Spearman rank correlation coefficients of the order of 
-0.1 to -0.3. This result shows that the [Oll]/Hj3 calibration 
given in Eq.[23]is not significantly sensitive to variations in the 
dust attenuation. This calibration can thus reliably been used in 
any sample without applying a correction. This conclusion is not 
true for the [Oll]/Hj3 calibration given in Eq. [24] which shows 
similar, but lower, residual slopes as compared to single-line cal- 
ibrations. 

Fig. [T6] shows the [Oll]/Hj3 calibrations applied on the two 
sub-samples defined in Sect. [2] with two different dust proper- 
ties: (t^^^) = 0.63 and _(Tp^i) = 1.84. It confirms that the 
calibration defined in Eq.[23]is reliable when applied on samples 
with different dust properties: the dispersion does not signifi- 
cantly changes, and the observed residual slopes are not signifi- 
cant (Spearman rank correlation coefficients of the order of 0.2). 
In both cases, the observed systematic shifts are in agreement 
with the correction formula given in Eq.[25]and the coefficients 
given in Table [T] 



4.3.2. Dependence on metallicity 

We also study the residuals of our new calibrations, as compared 
to the reference CLOl SFR, as a function of the gas-phase oxy- 
gen abundance estimated with the CLOl method. For the single- 
line calibrations. Fig. [T7] shows a significant correlation between 
these two quantities with Spearman rank correlation coefficients 
of the order of —0.5 to —0.6. Thus, we can derive linear rela- 
tions which may be used to correct these calibrations for a de- 



Table 2. Coefficients to be used to correct the SFR derived from 
an observed calibration, as a function of the gas-phase oxygen 
abundance following Eq. . (first three calibrations) or Eq. (two 
last calibrations). 



calibration 


a 


b 


c 


Ha 


-0.53 


4.80 




[On] 


-1.50 


13.35 




Hi8 


-0.74 


6.67 




[Oll]/Hj8 (Eq.[23l) 


1.47 


-25.18 


107.5 


[Oll]/Hj3 (Eq.[2l 


0.51 


-8.96 


39.3 



a=-0.037 b = -0.210 
spear=-0.144 
mean=-0.190 
-rms = 0.215 




a=-0.161 b 
spear=-0.318 
mean=0. 165 
= 0.279 



log SFR (Mg/yr) 



a=-0.080 b = -0.097 
spear=-0.290 
mean= -0.053 
-rms = 0.1B7 




log SFR (Mg/yr) 



a=-0.241 b=0.052 
spear=-0.530 
~^""n=0.006 
".263 




log SFR (Mg/yr) 



log SFR (Mg/yr) 



Figure 18. Same legend as in Fig. [2] The studied calibration is 
the new [Oll]/Hj3 calibration given by Eq. [23] (top) or Eq. [24] 
(bottom), applied on the two metallicity sub- samples defined in 
Sect.[2] Left: (12 + log(0/H)) = 8.68, right: (12 + log(0/H)) = 
9.07. 
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Figure 17. Same legend as in Fig.O except that the residuals of the SFR calibration are now plotted as a function of the gas-phase 
oxygen abundance estimated with the CLOl method. The studied calibrations are, from left to right: Ha, [Oil], Hj3, [Oll]/Hj3 
(Eq.[23]), and [Oll]/Hj3 (Eq.[24]). All calibrations are related to observed quantities. 



sired studied sample, for which one has a rough estimate of the 
mean metallicity. The corrective formula as derived from Fig.fTTl 
is: 

log (SFR"^"") = log (SFR) + a X jc + /? (26) 

with X = 12 + log(0/H). The a and b coefficients are summa- 
rized in Table [21 We note that this correction formula give an ap- 
proximately null correction for a metallicity equal to the mean 
value in the SDSS DR4 sample: 12 + log(0/H) = 8.89. 

For both [Oii]/Hj3 calibrations, we have fitted a second de- 
gree curve to the residual, which leads to the following corrective 
formula: 

log (SFR"^"") = log (SFR) -^axx^^bxx^c (27) 

with X = 12 + log(0/H). The a,b and c coefficients are summa- 
rized in Table [21 

For the [Oll]/Hj3 two-lines calibration given in Eq. [24l we 
see that the residuals do not show any more any strong signifi- 
cant correlation with the metallicity, with a Spearman rank cor- 
relation coefficients of 0.08. This result shows that the [Oii]/Hj3 
calibration given in Eq.[24lis not significantly sensitive to vari- 
ations in the metallicity. This conclusion is not true for the 
[Oll]/Hj3 calibration given in Eq. [211 which shows a Spearman 
correlation coefficient of 0.57. 

Fig. [T8l shows the [Oll]/Hj3 calibrations applied on the two 
sub-samples defined in Sect. [21 with two different metallic- 
ity properties: (12 + log(0/H)) = 8.68 and (12 + log(0/H)) = 
9.07. It confirms that the calibration defined in Eq. [211 is reli- 
able when applied on samples with different dust properties: the 
dispersion does not significantly changes. In both cases, the ob- 
served systematic shifts are in agreement with the correction for- 
mula given in Eq.[27land the coefficients given in Table [21 

We conclude that: 

- The [Oll]/Hj3 calibration given in Eq.[23lis weakly sensitive 
to variations in dust attenuation. It can be applied on sam- 
ples with an assumed metallicity, using the corrective for- 
mula given in Eq.[27l 

- The [Oll]/Hj3 calibration given in Eq.[24lis not sensitive to 
variations in metallicity. It can be applied on samples with an 
assumed dust attenuation, using the corrective formula given 
in Eq.[25l 

5. Conclusions 

We draw the following conclusions from our study: 

- As already shown in many previous studies, we confirm from 
SDSS DR4 data the the best emission-line calibration of the 



SFR is based on Ha\ this one being corrected for dust atten- 
uation (see Sect. 13.2.1b . This calibration has an uncertainty 
ofO.ndex. 

- If the dust has been estimated from the wrong assumption of 
a constant i ntrin sic B aimer ratio, the standard one-parameter 
[Kennicuti (119981) scaling law gives good enough results. 
Nevertheless, the SDSS DR4 data has shown a -0.11 dex 
systematic shift which can be corrected either by applying 
a mean shift in the opposite direction, or by using a two- 
parameter power law (see Fig. [21 Fig. [51 and Eq.[7l). 

- When Ha^ is not observed but a corr e ction for dust atten- 
uation is still available, the iKennicuti (1 19981) law ba sed on 
[Oii]^ (and corrected for dust by iKewley et al.ll2004l) gives 
good results but shows a higher dispersion (see Fig. [21). 

- This dispersion may be reduced by taking into account 
the dependence of [Oii]^ on the metallicity. Howeve r, the 
metallicity correction proposed by lKewley et al.l (|2004|) does 
not end to a significant improvement, while it relies on a 
very uncertain calibration of the metallicity, as shown by 
Kewley & ElHson (2008) (see Fig. [Hand Fig. [51). 

- A very good correction of this effect is rather obtained by 
using the [Oii]^ and Hj3^ lines together or, even better, Hj3^ 
alone (see Eq.[9l Eq.[Tniand Eq.[n]). 

- We caution the reader against the use of the inadequate SFR 
calibration when data is corrected for dust estimated from 
another method (e.g. CLOl method). If the method to derive 
dust attenuation is not biased towards the wrong assumption 
of a constant intrinsic B aimer ratio, then a different SFR cali- 
bration should be used (see Sect. [33H Fig.[6l Eq.[I2 Eq.[T3l 
and Eq.O. 

- We advise the reader not to use dust estimated from SED 
fitting in order to calculate SFR from emission lines. This 
method leads to highly uncertain results because of the high 
dispersion in the stellar-to-gas attenuation ratio (see Fig.[71). 
We emphasize that dust attenuation estimated from SED fit- 
ting would be reliable only if it is applied to the light emitted 
by stars, not to emission lines. 

- When no estimation of the dust attenuation is available, it 
is common in the literature to assume a mean correction, 
typically Ay = \. We advise the reader to use such method 
with great care for two reasons: (i) the SFRs recovered with 
the assumption of a constant dust attenuation are biased by 
a non-negligible residual slope, coming from the correla- 
tion between dust attenuation and SFR (see Sect. 14.1.11 and 
Fig. [SI); (ii) additional biases come from the choice of the 
assumed dust attenuation which is likely not to be the right 
one. The choice of a wrong assumed dust attenuation leads 
to non-negligible systematic shifts (see Sect. 14.1.21) . 
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- We have derived new direct calibrations between the SER 
and the observed line luminosities which are still quite poor 
in terms of dispersion, even if they do not show any more a 
significant residual slope (see Eig.O Eig.[TTJ Eq.[T71 Eq.fTSl 
andEq.[T9l). 

- The calibrations p roposed by iMoustakas et alJ (l2006l) or 
IWeiner et"aD ('2007) which include a correction based on the 
B-bmd or H-band ^-corrected absolute magnitudes do nei- 
ther significantly better nor significantly worse (see Eig.[T3]). 

- We derive two new two-lines calibrations based only on the 
[On] and Hj3 observed line flux. These calibrations give the 
best results among all calibrations based on observed quanti- 
ties (i.e. not reliably corrected for dust attenuation). It shows 
no systematic shift, no significant residual slope, and a re- 
duced dispersion (see Eig.O Eq.[23]and Eq.O. The mini- 
mum uncertainty with data not corrected for dust attenuation 
is 0.23 dex. 

- We have studied the relation between the residuals of our 
calibrations based on observed quantities, and the dust atten- 
uation (see Eig. [15]). There is a clear correlation for single- 
line calibrations, which can be corrected (see Eq. [25] and 
Table[T]). We have also found a correlation between the resid- 
uals of our calibrations based on observed quantities (see 
Eig. [17]), which can be corrected for single-line calibrations 
using Eq.[26]and Table[2] 

- Among our two new [Oii]/Hj3 calibrations, the one defined 
by Eq. [23] is designed to be used on a sample with an as- 
sumed metallicity and an unknown dust attenuation, while 
the one defined by Eq.[24]give better results on samples with 
unknown metallicity and an assumed dust attenuation. 
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